BU-CCS-960601 



Interface Roughening in a Hydrodynamic Lattice-Gas Model 

with Surfactant 

Francis W. Starr, Stephen T. Harrington, Bruce M. Boghosian, and H. Eugene Stanley 
Center for Polymer Studies, Center for Computational Science, and Department of Physics, Boston University, Boston, MA 

02215 USA 
(LT5906: Received 18 June 1996) 



0\ 
ON 
t— I 

Oh! 

<D 

o 

(N 

> 
00 

oo 

o 



i 

c 

o 
o 



13 



Using a hydrodynamic lattice-gas model, we study interface growth in a binary fluid with various 
concentrations of surfactant. We find that the interface is smoothed by small concentrations of 
surfactant, while microemulsion droplets form for large surfactant concentrations. To assist in 
determining the stability limits of the interface, we calculate the change in the roughness and 
growth exponents a and (3 as a function of surfactant concentration along the interface. 

PACS numbers: 68.10.-m, 05.50.+q, 47.11. +j 



The study of rough interfaces is of important exper- 
imental and theoretical interest and has wide interdis- 
ciplinary applications, including fluid imbibition experi- 
ments, flow in porous media and fluid-fluid displacement 

. Rough interfaces have been extensively studied by di- 
rect integration of continuum equations and using simple 
discrete lattice models. The development of lattice-gas 
models makes it possible to study interface roughening in 
hydrodynamic systems. Lattice-gas models have evolved 
from simple one-component Navier-Stokes fluids to 
two-component immiscible lattice-gas (ILG) models jj), 
and, most recently, to a model including amphiphilic par- 
ticles [f|. Lattice-gas models can reproduce the dynam- 
ics on mesoscopic scales, allowing for the investigation of 
non-equilibrium behavior over a much broader range of 
length and time scales than is possible with molecular 
dynamics 0. The simplicity of the collision rules, the 
exact conservation of mass and momentum, and the nat- 
ural underlying kinetic fluctuations suggest that lattice- 
gas models are an appropriate choice for studying the 
scaling behavior of interfaces in complex hydrodynamic 
systems. 

In this Letter, we study interfacial roughening in the 
presence of a surfactant by including amphiphilic parti- 
cles that tend to reside at the binary fluid interface ||. 
Flekk0y and Rothman recently studied fluctuating inter- 
faces in an ILG model without surfactant and described 
the scaling properties of these interfaces 0. When sur- 
factant is added to an initially flat interface, we find that 
the roughness exponent a and growth exponent [3 de- 
crease for small concentrations of surfactant |^|; however, 
for larger concentrations of surfactant, (3 increases while 
a continues to decrease. The continued addition of sur- 
factant causes the interface to spontaneously break up 
into a microemulsion phase. We also find that the sat- 
urated width of the interface has a minimum when the 
surfactant concentration along the interface matches the 



binary fluid concentration in the bulk. 

Far from any interface, the present model approaches 
the one-component FHP-II model m , which employs a 
triangular lattice with each lattice site occupied by up 
to seven particles, each of which is either at rest or has 
unit velocity directed along one of the six lattice vectors 
toward neighboring sites. No two particles may have the 
same velocity at a lattice site. The rules dictate that a 
particle propagates along a lattice direction until a col- 
lision occurs. When no surfactant is present, the cur- 
rent model reduces to the two-component ILG model || , 
which employs two types of particle and assigns each of 
the two species a color, "red" or "blue". Each collision 
must conserve red particles, blue particles, and the total 
momentum. Additional collision rules are included that 
give rise to the aggregate behavior of an immiscible fluid. 
These rules can be described using an electrostatic anal- 
ogy |Q, where the color flux J(x, t) of an outgoing state 
is the difference between the red and blue momenta. The 
color field E(x, t) is the gradient of color between neigh- 
boring sites. At each time step, the color work - defined 
by the inner product J • E - is minimized at each lattice 
site, which has the effect of inducing phase separation. 

With the addition of surfactant molecules, the ILG be- 
comes a three-component model having the potential to 
model amphiphilic systems To make the behavior of 
the surfactant consistent with a molecule composed of 
hydrophobic and hydrophilic ends, the surfactant is rep- 
resented as a "color dipole" with a dipole vector <j(x,i) 
and produces a dipolar field P(x,i) at neighboring sites. 
A term proportional to a ■ E, representing the interaction 
of colored particles with dipoles, is added to the Hamil- 
tonian. Similarly, we add a term proportional to a ■ P 
to include the dipole-dipole interaction. Minimization of 
the modified Hamiltonian favors the surfactant particles 
lining up along the interface between red and blue parti- 
cles and makes it unfavorable for surfactant particles 
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FIG. 1. Typical configurations (a) for the L = 32 case 
without surfactant (p s = 0); the concentration in each region 
is p — 0.5; (b) for a surfactant concentration p s — 1.0 along 
the interface; binary fluid particles are replaced by surfactant 
particles so that the total particle number remains unchanged; 
(c) for the micellular microemulsion phase (p a — 1.5). Note 
that the periodic boundary conditions require two interfaces. 



to align with each other at neighboring sites, modeling 
the expected behavior of an amphiphilic molecule |10|. 

We perform simulations of the three-component model 
using a sequence of two-dimensional lattices of edge 
L = 16, 32, 64, 128 and 256 with periodic boundary con- 
ditions in both directions. Initially, red and blue particles 
are placed on the lattice separated by two flat interfaces 
which are evenly spaced and oriented perpendicular to 



one of the lattice directions (T^| , as shown in Fig. 
In each region the particle concentration is p = 0.5 [121. 
Each particle, with uniform probability, is either given 
zero velocity or unit velocity in one of the six lattice di- 
rections. Surfactant molecules are placed between the 
interfaces with varying concentration and angular orien- 
tation randomly selected between between and 27r. (Af- 
ter a single collision step, the angular orientation of the 
surfactant molecules tends to align perpendicular to the 
interface, as this is energetically favorable.) The system 
then evolves according to the collision rules of the model. 
Our interest focuses on how the addition of surfactant to 
the interface mediates the roughening process. 

We first calculate the average interface width 
W(L,t) = (/i 2 (x,<) - /i(x,i) 2 } 1/2 at logarithmically 
spaced time intervals with no surfactant present. We find 
that W(L,t) ~ t 13 until a crossover time t x (L), where- 
upon it saturates at a value W sat {L) ~ L a [0,0, with 
P = 0.33 ± 0.02 and a = 0.50 ± 0.02 (Fig. |). Hence z = 
a/ (3 = 1.50 ± 0.03, where z is the dynamical exponent, 
defined by t x (L) ~ L z . In Fig. ||(b), we plot the rescaled 
width W(L,t)/L a against the rescaled time t/L z and 
confirm the scaling hypothesis, W(L,t) ~ L a f(t/t x ) 
0. Our results match the scaling exponents a = 1/2 
and (3=1/3 for the Kardar-Parisi-Zhang (KPZ) uni- 
versality class |l3| , and are comparable to earlier results 
obtained by Flekk0y and Rothman . 

Next we consider interface roughening with surfactant 
present. In the three-component model, for surfactant 
concentrations below the critical concentration for spon- 
taneous emulsification, the surfactant particles tend to 
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FIG. 2. (a) Interface growth and roughening of a bi- 
nary fluid for system sizes L = 16(o), 32(D), 64(0), 128(A) 
and 256(<l). At early times, the width W ~ t 13 where 
/3 = 0.33 ± 0.02. (b) Rescaling to test the scaling hypoth- 
esis and the values of a and f3. Error bars are approximately 
equal to the symbol size. The inset shows the size scaling of 
saturation width, where a = 0.50 ± 0.02, in agreement with 
the values predicted in ref. pj. 



line up along the interface between the two fluids. Since 
we are interested in the role that the surfactant plays in 
the kinetic roughening process, we start with an initial 
configuration where the fluids are separated and we re- 
place red and blue particles with surfactant particles in 
a single row along the interface (Fig. |l](b)). We define 
the surfactant concentration p s = N s /L, where N s is the 
number of surfactant particles placed along the interface. 
To obtain p s > p, we replace red and blue particles with 
surfactant particles on multiple rows [flif . 

In Fig. ||(a), we plot W(L,t) for all system sizes and 
p = p s = 0.5. We still see a growth regime that crosses 
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FIG. 3. (a) Interface growth and roughening for the three 
component model with a surfactant concentration p a = 0.5 
along the interface. System sizes and symbols are the same 
as Fig. ^ For t < t x , the growth rate is slower than for 
the binary case, and the saturation width is smaller. We find 
/3 = 0.27 ± 0.02. (b) Rescaled plot demonstrating data col- 
lapse. Error bars are approximately equal to the symbol size. 
The inset shows the roughness exponent a = 0.43 ± 0.02. 

over to a saturated regime. Compared to Fig. the 
case without surfactant, the growth rate of the width is 
slower and W sat is smaller. We find (3 = 0.27 ± 0.02 
and a — 0.43 ± 0.02. Fig. ||(b) confirms the scaling hy- 
pothesis using z — a/ (3 — 1.59 ± 0.03. These values for 
a and (3 are not the accepted values for the roughness 
and growth exponents of the KPZ equation in 1 + 1 di- 
mensions. The KPZ exponents also obey the sum rule, 
z + a — 2, a consequence of Galilean invariance jy. We 
obtain z + a = 2.02 ± 0.04 for this case. This result is 
surprising since lattice-gas models are not Galilean in- 
variant, due to the fact that the lattice constitutes a pre- 
ferred Galilean reference frame M. On the other hand, 




N s /L 

FIG. 4. (a) Dependence of the growth exponent /3 on sur- 
factant concentration p s = N s /L for a 32 x 32 system with 
p = 0.5. For small surfactant concentrations, the growth 
rate is reduced. j3 increases for p s > p until spontaneous 
emulsification occurs, (b) Dependence of saturation width on 
surfactant concentration for p = 0.5. For p s < p the surfac- 
tant reduces the value of the saturated width. For p 3 > p the 
interface becomes rougher as it approaches spontaneous emul- 
sification. (c) Wsat dependence on p 3 with a binary fluid con- 
centration p = 0.7. We again see that the behavior changes 
at p a = p. 



the flow equations of the lattice-gas model do contain a 
non-linear inertial term, as does the KPZ equation. 

We next consider the effect of the surfactant concen- 
tration on the growth of the interface and the saturated 
width. In Fig. [|(a), we show the dependence of f3 on p s 
for a lattice of size 32 x 32. We observe a dramatic re- 
duction in the growth exponent upon the addition of any 
surfactant to the interface. Additional surfactant causes 
little change in j3 up to p s = p, where (3 = 0.22±0.02. For 
p s > p, the growth exponent increases but never returns 
to the value of the pure case. For p s > 1.5, the inter- 
face spontaneously breaks up, giving rise to a miccllular 
microemulsion phase where the width is no longer well 
defined (Fig. 1(c)). The effect of p s on W sa t is similar 
to its effect on (3. For p s < p we see that the interface 
width saturates at a value that decreases for increasing 
p s (Fig. ||(b)). However, for p s > p, W sa t increases until 
we reach the microemulsion phase. 

We can understand the behavior of W sa t by considering 
the effect the surfactant has on the binary fluid interface. 
From the Hamiltonian, we know that the color-dipole en- 
ergy, (T-E, required for surfactant to be removed from the 
interface is large compared to the energy for the surfac- 
tant to remain on the interface. As a result, the surfac- 
tant particles are effectively "anchored" at the interface. 
For p s < p, the surfactant reduces the fluctuations of 
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the interface because it is anchored in position between 
red and blue particles. As a result, W sa t is smaller. For 
p s > p, the color-dipole interaction still forces surfactant 
particles to remain on the interface, but now surfactant 
particles are more likely to neighbor each other. Since the 
dipole-dipole interaction, a-P, makes it energetically un- 
favorable for dipole molecules to align with other dipole 
molecules at neighboring sites, the surfactant creates ad- 
ditional surface by roughening the interface to reduce the 
number of neighboring surfactant molecules, thereby in- 
creasing W sa t- To demonstrate that the turning point 
in W S at depends on our choice of p, we measure W sa t 
for various p s with a binary fluid concentration p = 0.7 
(Fig. H(c)). We again see W sa t decreasing for p s < p and 
increasing for p s > p. 

To investigate the effect of p s on a, we consider sys- 
tem sizes L = 16, 32, 64 and 128 for p s = 0.8 and 1.0. At 
p s = 0.8 we find f3 = 0.27±0.02 and a = 0.36±0.02, while 
for p s = 1.0 we obtain (3 = 0.30±0.02 and a = 0.23±0.02. 
Using the exponents for these concentrations, we do not 
find z ± a = 2. In Table I we summarize the behavior 
of the exponents a, (3, z, and z + a as a function of the 
surfactant concentration. We notice that a decreases as 
we increase the surfactant concentration. 

In summary, we have found that a three-component 
lattice gas can be used to test the effect of surfactant of 
the scaling properties of hydrodynamic interfacial growth 
and roughening. In particular, we have found that sur- 
factant alters the growth and roughness exponents, with 
a decreasing and increasing as we increase p s above 
the concentration of the binary fluid. These results may 
assist in understanding the stability limits of interfacial 
growth in the presence of a surfactant and the conditions 
for the formation and breakup of microemulsion droplets. 



TABLE I. Roughness, growth and dynamic exponents 
for various surfactant concentrations measured from the 
simulations. Also included are the number of samples that 
have been averaged. 



0.0 
0.5 
0.8 
1.0 







0.50 ±0.02 0.33 ±0.02 1.50 ±0.03 2.00 ± 0.04 

0.43 ±0.02 0.27 ±0.02 1.59 ±0.03 2.02 ± 0.04 

0.36 ±0.02 0.27 ±0.02 1.33 ±0.03 1.69 ±0.04 

0.23 ±0.02 0.30 ±0.02 0.77 ± 0.03 1.00 ±0.04 



Samples 



170 
170 
150 
150 
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